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ABSTRACT 

Using the Teukolsky and Sasaki-Nakamura equations for the gravitational per- 
turbation of the Kerr spacetime, we calculate the post-Newtonian expansion of 

•rH . 

X ■ the energy and angular momentum luminosities of gravitational waves from a test 

particle orbiting around a rotating black hole up through P 5 / 2 N order beyond the 
quadrupole formula. We apply a method recently developed by Sasaki to the case 
of a rotating black hole. We take into account a small inclination of the orbital 
plane to the lowest order of the Carter constant. The result to P 3 / 2 N order is 
in agreement with a similar calculation by Poisson as well as with the standard 
post-Newtonian calculation by Kidder et al. Using our result, we calculate the 
integrated phase of gravitational waves from a neutron star-neutron star binary 
and a black hole-neutron star binary during their inspiral stage. We find that, 
in both cases, spin-dependent terms in the P 2 N and P 5 / 2 N corrections are im- 
portant to construct effective template waveforms which will be used for future 
laser-interferometric gravitational wave detectors. 



1. Introduction 



The last stage of an inspiraling compact binary such as binary neutron stars 
is one of the promising sources of gravitational waves for the near-future laser- 
interferometric detectors such as LIGO [11 and VIRGO.' 21 The main reasons why it is 
most promising are that the event rate is expected to be ~ 3 events/yr within 200 
Mpc from a statistical study and that it is possible to theoretically predict the 
amplitude of gravitational waves with good accuracy. Such a binary is, moreover, 
not only a strong source of gravitational waves, but also has a possibility to be- 
come a treasury of physics of neutron stars) 41 cosmology, [51 theories of gravity) 61 etc., 
provided we obtain the data about binaries such as masses, spins, distance to the 
earth, and so on. However, to obtain those data with sufficient accuracy, it is nec- 
essary to construct theoretical template waveforms whose phasing has a fractional 
accuracy of less than 10~ 4 .' 41 Hence considerable efforts have been paid to construct 
such theoretical templates! 7-151 

To calculate gravitational waves from an inspiraling compact binary, the stan- 
dard method employed is the post-Newtonian expansion (PNE) of the Einstein 

ryi 

equations, in which the equations are expanded in terms of a small parameter 
v ~ (M/r) 1 / 2 , where M and r are the total mass and orbital length-scale of the 
system, respectively. Despite much effort, however, calculations have been success- 
ful to only a few orders in v beyond the leading (Newtonian) order so far. More 
fundamentally, the nature of PNE has not been clarified due to its complexity; 
nobody knows the convergence property of PNE nor the validity of the polyno- 
mial expansion in t>. [161 Given this situation, it is highly desirable to have a method 
which is complementary to the standard PNE. The perturbative study of a black 
hole spacetime is one of such, in which we consider gravitational waves radiated 
by a particle of mass fi <C M orbiting around a black hole. This method, though 
restricted to the case of /i <C M, is very powerful because we can calculate fully 
general relativistic corrections of gravitational waves by means of relatively simple 
analyses. It is then fairly straightforward to evaluate the post-Newtonian correc- 
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tions which are to be calculated in the standard PNE. This direction of research 
was first done analytically by Poisson [81 to P 3 / 2 N order and numerically by Cutler 
et al. 191 to P 3 N order. Then a highly accurate numerical calculation was done by 
Tagoshi and Nakamura' 141 to P 4 N order and an analytical calculation to the same 
order was done by Tagoshi and Sasaki' 151 which confirmed the result of Tagoshi and 
Nakamura. In particular, the appearance of logf terms in the energy flux at 0(v 6 ) 
(P 3 N order) and at 0(v 8 ) (P 4 N order) is confirmed and it is clarified that the 
accuracy of the energy flux to at least P 3 N oder is necessary to construct template 
waveforms for optimal use of the interferometric data. 

In this paper, we extend the analysis of Tagoshi and Sasaki ' 1j1 to the case 
of circular orbits around a rotating black hole to see the effect of spin. Some 
analytical' 10 ' 121 and numerical analyses' 11 ' 131 of its effect at its leading order (i.e., to 
P 3 / 2 N order) have been performed. They found that it gives rise to a large error in 
the integrated phase of gravitational waves from inspiraling binaries for a typical 
value of the spin angular momentum if it is ignored in the template. The large 
effect of the P 3 / 2 N order spin-orbit term causes us great anxiety about the effects 
of higher order terms due to the spin angular momentum since the convergence 
property of the PNE has been found to be slow.' 4 ' 14 ' 10 ' Hence we study the effect of 
spin to P 5 / 2 N order in this paper. 

The paper is organized as follows. In section 2, we show the basic formulas 
to perform the PNE in our perturbative approach. First we show the PNE of 

[171 T181 

the Teukolsky radial function using the Sasaki-Nakamura equation. We must 
also perform the PNE of the angular equation, which is given in Appendix D. 
Then we describe the PNE of the source terms. We consider circular orbits around 
a Kerr black hole, that is, those at constant Boyer-Lindquist radial coordinate. 
These circular orbits are, however, not necessarily on the equatorial plane and the 
emission rates of gravitational waves are different for different orbital inclination 
angles (or more precisely, different values of the Carter constant).' 131 To see the 
leading effect of the orbital inclination, we consider orbits with small inclination 
angles. In section 3, the energy and angular momentum luminosities to 0(v 5 ) 
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beyond Newtonian are derived. In section 4, using the results given in section 
3, we estimate errors in the accumulated phase of gravitational waves caused by 
the non-vanishing spin angular momentum. Section 5 is devoted to the summary. 
Throughout this paper we use the units of c = G — 1 and GMq/c 2 = 1.477km. 

2. Formulation 

To calculate gravitational radiation from a particle orbiting around a Kerr black 

[17] 

hole, we start with the Teukolsky equation. We focus on the radiation going out 
to infinity described by the fourth Newman- Penrose quantity, ^4 j 191 which may be 
expressed as 

/ e irrvfi 
due-"* 77^ -2SZ(8)RlmM, (2.1) 

£,m 

where -2<S^w is the spheroidal harmonic function of spin weight s = —2, which are 
normalized as 



7T 



2 Sr| 2 sin0d0 = l, (2.2) 



o 



and A is the eigenvalue. The radial function Ri muJ {r-) obeys the Teukolsky equation 
with spin weight s = — 2, 

A2 ^^r") - V(r)Re ™ = T ^ (r) ' (2 - 3) 

where Tg mu} (r) is the source term whose explicit form will be shown later, and 
A = r 2 - 2Mr + a 2 . The potential V(r) is given by 

K 2 + Aiir — M)K n . 
V(r) = y — '— + 8iur + A, (2.4) 

where K = (r 2 + a 2 )u — ma. 
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The solution of the Teukolsky equation at infinity (r — > oo) is expressed as 



00 



^ mw(r) m£ y rfr — av) — = ^ r e ' (2 - 5) 



where r_|_ = M + V M 2 — a 2 denotes the radius of the event horizon and -Rf^ is 
the homogeneous solution which satisfies the ingoing-wave boundary condition at 
horizon, 

where k = u — ma/2Mr + and r* is the tortoise coordinate defined by 



dr* r 2 + a 2 
~dr ~ A ' 



(2.7) 



For definiteness, we fix the integration constant such that r* is given explicitly by 

dr* , 
— — or 

(2 8") 

2Mr + , r - r + 2Mr_ , r - r_ ' 1 ' ; 

= r H — In In 

r + - /• 2M r + - /• 2.1/ 

where r± = M ± VM 2 - a 2 . 

Thus in order to calculate gravitational waves emitted to infinity from a particle 
in circular orbits, we need to know the explicit form of the source term T^ mw (r), 
which has support only at r = ro where ro is the orbital radius in the Boyer- 
Lindquist coordinate, the ingoing-wave Teukolsky function -R^i W (r) at r = ro, 
and its incident amplitude Bj 7 ^ at infinity. We consider the expansion of these 
quantities in terms of a small parameter, v 2 = M/r-Q. Note that v is approximately 
equal to the orbital velocity, but not strictly equal to it in the case of a ^ 0. 
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In addition to these, we need to expand the spheroidal harmonics and their 
eigenvalues in powers of auj. Since uj = 0(Q) where Q is the orbital angular velocity 
of the particle, we have auj = O(Muj) = 0(t> 3 ). Thus the expansion in powers of 
auj is also a part of PNE. Note also that the spin parameter of the black hole a 
does not have to be small but can be of order M. We defer this expansion of the 
spheroidal harmonics to Appendix D. 

In the following, we consider the PNE of the ingoing-wave Teukolsky function 
and the source term separately. Since the angular eigenvalues A comes into play in 
the radial equation, we set 

\ = \o + auj\ l + a 2 uj 2 \ 2 + 0(v % ), (2.9) 

where A = {£ - l){£ + 2) and Ai = -2m(£ 2 + £ + A)/(f + £) have been obtained 
previously) 201 and A2 is derived in Appendix D. 

2.1. Homogeneous solution 

The PNE of the Teukolsky equation for a rotating black hole was first per- 
formed by Poisson up to 0(t> 3 ) beyond the quadrupole formula! 121 In his method, 
he introduced a parameter e = 2Mco ~ 0(t> 3 ) and expanded the equation in terms 
of it to the first order in e, since only the homogeneous solution to 0(e) is needed to 
obtain the energy flux up to 0(v 3 ) level. The purpose here is to obtain the higher 
terms, 0{v n ) (n > 3). Specifically in this paper we consider the gravitational wave 
luminosity to 0(v 5 ). Hence we need the homogeneous solution accurate up to 
0(e 2 ). However, it is very difficult to treat the Teukolsky equation itself to obtain 
the higher order corrections because the Oth order solutions (i.e., the kernel func- 
tions) already have a quite complicated form! 121 Hence instead of the Teukolsky 

[18] 

equation, we use the Sasaki-Nakamura (SN) equation, which is obtained by a 
certain transformation of the Teukolsky equation. The reason to use it is that the 
SN equation is a generalization of the Regge- Wheeler (RW) equation for a = 
to the case of a 7^ and hence has a more tractable structure. In addition, we 
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can make use of algebraic formulas for the PNE of the RW equation developed by 
Poisson [81 and by Sasaki! 211 

The SN equation has the form, 



(2.10) 



where the explicit forms of F and U are given in Appendix A. It is obtained by 
the transformation from R^ mL0 to Xg mijj as 



it -\u I x 



:2.1V. 



where Ximcu — Xi mu; A/(r 2 + a 2 ) 1 / 2 , and the functions a, (3 and rj are shown in 
Appendix A. Conversely, we can express Xi muJ in terms of Ri muJ as 



X tmu = (r 2 + a 2 ) 1 / 2 r 2 J-J- 



(2.12) 



where J_ = (d/) — i(K/A). Then the ingoing- wave solution X^^ which corre- 
sponds to nas the following asymptotic behavior, 



a out „iu>r* , \in a —iiar* 
/~t —ikr* . 



(2.13) 



r — > — oo, 



where Aff^, A^ and Q mw are respectively related to Bf^, and D 

defined in Eq.(2.6) as 



'lmu> 



r>in 
n lmuj 



r>out 
n £muj 



A in 

4a; 2 



4u< 
1 



A OUt 



(2.14) 



7 



where Co is given in Eq.(A.3) of Appendix A and 

d emi0 = ^2Mr+[(8 - 2AiMuj - 16M 2 u 2 )r 2 + 

+ {Yliam - 16 M + IQamMcu + 2AiM 2 cu)r + - Aa 2 m 2 - YliamM + 8M 2 ]. 

(2.15) 

Now let us consider the PNE of XI™.. It is characterized by the ingoing- wave 
boundary condition at horizon, X^^ ~ e~ ikr as r* — > — oo. However, since r* 
cannot be expanded in terms of e at r* — > — oo, a naive expansion of the SN 
equation in terms of e would obscure the boundary condition at horizon. One 

[21] 

prescription to circumvent this difficulty has been suggested by Sasaki in the 
case of the Schwarzschild black hole, namely to separate out the factor e~^ r 
from Xff^ from the beginning. Here we generalize this method to the case of the 
Kerr black hole. 

First we introduce the variable z = ur and 



* i 
z = z + e 



\n(z — z+) \n(z — Z-) 



Z + — Z- Z + — Z- 

= ur* + e In e , 



(2.16) 

where z± = ur±. For later convenience, we also introduce a non-dimensional 



parameter q = a/M. Hence, for example, au = \qe. z± = \e(l ± a/1 — q 2 ). As 
mentioned before, q is not necessarily very small but can be of order unity. 

Next we define a function (f>(z) as 

4>(z) = I dr ( — — uj 



' vA 7 , (2.17) 

* e 1 , z — z + 
= z — z mq In , 

2 z + — Z- z — z- 

which generalizes the phase function uo(r* — r) of the Schwarzschild case. Note that 
e -#(«) e -ikr* ag r * _^ _ QO an( j ^ e -tuj(r*-r) ag r * _^ +QO> Then we set 

XZ U = V z 2 + a 2 ujH em (z) exp (-i<f>(z)) . (2.18) 
By this prescription, it is easily to implement the ingoing-wave boundary condition 
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Inserting Eq.(2.18) into Eq.(2.10) and expanding it in terms of e = 2Mu, we 
obtain 

L(°)[£ £m ] = eL^ em ] + eQ«fe m ] + e 2 Q^[^ em ] + 0(e 3 ), (2.19) 
where L^°\ Q^> and are differential operators given by 

_ 1 d 2 + / 1 + 2i \ d _ / 4 i + 1\ 

zdz 2 yz 2 z J dz \z 3 z 2 z J 
n (D _ »gAi d 4img 4mg Aig + 2mg 

^ 2^ ^ /(/ + 1)^3 + /(/ + l)z 2 + 2 Z 2 ' 1 J 

for arbitrary £ and 

( 2 ) <? 2 d 2 (9q 2 - 2m 2 q 2 - 6imq i(2m 2 q 2 - 9q 2 ) - 18mq\ d 
^ =_ 4i2^2~V 36? + 108^2 ) d~z 

( 9q 2 + 5m 2 q 2 - 3imq i(2m 2 q 2 - 9q 2 ) 
+ V 9? + 54? 

-Q3imq + 81q 2 - 2m 2 g 2 
378^2 

(2.23) 

for £ = 2. We shall see that we do not need for I > 3 for the post-Newtonian 
order we consider in this paper (i.e., up to 0(v 5 ) beyond Newtonian). Expanding 
£f m in terms of e as 

oo 

= (2-24) 

n=0 

we obtain from Eq.(2.19) the iterative equations, 

L<»)[^]=0, (2.25) 

£ <0) K£] = £ <I) K2]+e (1) K2]-< ) , (2.26) 

i (0, [^l = ^'KS] + 0">[«ai + (2) [^] - <>■ (2-27) 
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The general solution to Eq.(2.25) is immediately obtained as 

t (o) _ (0) , +fl (0) 



(2.28) 



Now we consider the boundary condition. The argument is parallel to that 
given in Ref.[21] for the Schwarzschild case. The condition that Xff^ ~ e -ikr ag 
r* — > — oo implies \J z 2 + e 2 q 2 /4£ i g m (z) should be regular at z — oor+ ~ e. Since 
e can be made arbitrarily small, \/ z 2 + e 2 q 2 /A^^(z) should be no more singular 
than 0(z~ n ) at z — 0. Thus in particular we have = af^ ji- For convenience, 
we set Oip = 1. Taking into account of the behavior of the lowest order solution, 
we then infer that \/ z 2 + e 2 q 2 /A^j i (z) must be no more singular than z e+1 ~ n at 

(n) 

z = 0. Hence for n < 2, the boundary condition is that Q m ; is regular at z — 0. 

As noted previously, in the case of a circular orbit of radius tq, the source term 
Timu nas support only at r = tq and uotq = 0(Qro) = 0(v). Hence we only need 
If at z — 0(v ) <S 1 to evaluate the source integral, apart from the value of the 
incident amplitude A^^. Hence the PNE of X 1 ^^ corresponds to the expansion 
not only in terms of e = 0(t> 3 ) but also z by assuming e < z < 1. In order to 
evaluate the gravitational wave luminosity to 0(v 5 ) beyond the leading order, we 
must calculate the series expansion of £^ in powers of z for n = to £ = 4, for 
n = 1 to i = 3 and for n = 2 to I = 2 (see Appendix C). On the other hand, the 
accuracy of A^^ we need for this purpose is O(e) or n < 1. Thus we need the 
asymptotic behavior of £^ at infinity only for n = 1 (that for n = is trivially 
obtained). 

To calculate ii mu] to the accuracy discussed above, we rewrite Eqs.(2.26) and 
(2.27) in the indefinite integral form by using the spherical Bessel functions jt an d 

ni, 

z z 

dm =ntj dzz 2 Je W^ - it J dzz 2 n e W^ (n = 1, 2). (2.29) 
The series expansion formulas for £^ around z = is easily obtained if we know 

in) 

those for W)'. Then it is straightforward to impose the boundary condition at 
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z = 0. Further, if the integrations can be done exactly in closed analytic form, it 
is easy to extract out the values of the incident amplitude ^H" by examining the 
asymptotic behavior at infinity of thus obtained functions. As we shall immediately 
see, this is easily done for n = 1, which is sufficient for our purpose. 

For n — 1, if we set q — 0, Eq.(2.26) becomes the same equation as that 
discussed in Ref.[21]. Hence the only correction is to include the contribution from 

t(°)i 

to yield 



the term. Using the formulas developed in that paper, this is easily done 



(2.30) 



(i) = (i) (£-!)(£ + 3) . l 2 -4 . 

a e Je+ 2(£+l)(2£+l) 3e+1 2£(2£+l) 3e ~ 1 

+ ng{Q\2z — 7 — In 2z) — ji$\2z + ijt In z 

imq ( + 4 \ . imq ( (£+l) 2 + 4 \ . 
+ ~Y \i 2 (2£+l) ) Jl ~ l + ~Y V(^ + l) 2 (2£ + l) ) 31+11 

where C\{x) = — f x dtcost/t and Si(x) = f Q dtsmt/t are cosine and sine inte- 
gral functions, 7 is the Euler constant, and Oip is an integration constant which 
represents the arbitrariness of the normalization of Xff^. We set = for 
simplicity. 

In the actual calculation of the gravitational radiation to infinity, we need to 
know X l £ UJ {z) at e <C z <C 1. Using Eq.(2.17), we obtain the expansion of Eq.(2.18) 



as 



x tL) = V z 2 + a 2 u 2 & TO exp (-iip) 



+ e 2 



(2.31) 



. if (°) _ U {1) ma + V 0) (--ma + £ - - 
% ^lm 2^ £m q Z em \ 4 8 8 / 2 lm 



For n = 2, expanding Eq.(2.30) in terms of z and inserting it into Eq.(2.29), 
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(2) 

we have the series expansion of Q^. Inserting it into Eq.(2.31), we obtain 

in _ z 3 z 5 z 7 9 /imqz 2 13z 4 llimqz 4 6 \ 

- 15 - 210 + 7560 + ° (Z ) + f l^0~" W ~ ^780~ + ° ( ~ 'J 

+ e 2 C 2 + 2 7 2 7 mV .-^ + 0(^)), (2.32) 

- 4 - iso + 0(/) + 'Hs + IP + °( 25 ))' < 2 - 33 ' 

= ^ + 0(^ 7 ). (2.34) 

Once we have the homogeneous solutions of the SN equation, we have only to 
perform the transformation (2.11) to obtain the corresponding solutions of the 
Teukolsky equation. The results are 

j-fin z i 5 Hz® i j 23 z^ i g , , 
Wi * 2mw = 30 + 45 2f " 1260 ~ 420 Z + 45360 + 11340 2 

f—z 3 i o >? a mqz 4 

+ e mqz z H 

V 15 60 y 60 45 



15 60 * 60 45 
Alz 5 277i = 31 i fi 7mqz 6 \ 

— ~\~ 777. (7 J6 — £ — ) 

3780 22680 y 3780 1620 / 

2 /z 2 i o q 2 z 2 m 2 q 2 z 2 
+ £ v30 + 40 mgZ + "60 24b- 



i 



z 3 - + ± q 2 z 3 -J- m 2 q 2 z 3 ) , (2.36) 

30 90 120 / v ; 



60 

= ^1 + - -i! 

3mw 630 1260 3780 16200 

+ <252 - 189b m ^ - 756 " + ^2680-)' ^ 

z 6 iz 7 

URT ™ = 11340 + 28350' (2 ' 38) 



Here, it is worth noting that the terms linear in q{= a/M) at 0{y 2 ) beyond the 
leading term in each -Rl" are pure imaginary. This implies there is no linear term 
in q at the P 1 N order of the luminosity. Such terms will appear at the P 3 / 2 N order, 

[12] 

which is in fact what was found by Poisson. Further, it is expected that terms 
linear in q will not appear at the P 2 N order, but at the P 5 / 2 N order. By the same 
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argument, terms quadratic in q are expected to appear at the P 2 N order, but not 
at the P 5 / 2 N order. 

Next, we consider A^j^ to 0(e) — O(f^). The procedure is the same as that in 
the Schwarzschild case.' 211 Using the relations ~ —je-i ~ (—^Y +rln 2n-£, etc., 
we obtain the asymptotic behavior of ^ n at z = oo as 

~ -\h + - 111 Z ) n t + ^ 111 *» ( 2 - 39 ) 

where 



£(£+1) 

and ■?/>(£) is the digamma function, 

£-1 



fc 
fc=i 



To consider the asymptotic form of X t £^ w , we set £^ = fj^ + ijglaz and 
express the asymptotic form of fjp at z — > oo as 



f (i) _. p (i), , 



where h^jp and h!p are the spherical Hankel functions of the first and second kinds, 
respectively, which are given by 

hf\z) = j,{z)+in,{z) -> (-i) e+1 — , fc< 2 > (z) = j^)-m^) -> — . (2.41) 

^ ,2- 



From Eq.(2.39), we have 
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Then the asymptotic form of Xff^ becomes 



~ z 

1 

~ — e 
2 



h + <f^+ijelnz) + 



cxp 



—i\z — z — 



2y/l 



mq z — z-i 
m ■ 



7 2 z - Z- 



l + e(pW + i(Q)£ + \nz)) 
(!) 



1 



From this equation, -A^^ can be easily extracted out: 

1 



4 in _ %-£+l — ielne 



l + e(-| +i9 £) + 



(2.43) 



(2.44) 



where we use the fact that our* = z* — elne from our definition of z*. Specifically 
for I = 2, 3 and 4, to the orders respectively required, we have 



l 2mu) 



7T 



= -- 1 - e- + it [ - - 7 - log 2 + -le + 0(e 2 ), (2.45) 



4 m 



2 V3 

If, 7T . /13 



) + ^} + o (e2 ), 



(2.46) 
(2.47) 



Finally, the corresponding incident ampliutudes for the Teukolsky function 

are obtained from Eq.(2.14). 

2.2. The source term 



A test particle obeys the equations of motion, 



d9 

E— 

dr 

dt 
dr 



C — cos 2 6\ a 2 (l — E 2 ) 



It 



sin 1 6 



1/2 



= e(fl), 



-(aE 





sin 2 


9) 


lz 




sin 2 





r2 + a2 (E(r 2 + a 2 )-al z ), 



(2.48) 



where E, l z and C are the energy, the z-component of the angular momentum and 
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the Carter constant of a test particle, respectively, £ = r 2 + a 2 cos 2 6 and 

R = [E{r 2 + a 2 ) - al z } 2 - A[(Ea - l z ) 2 + r 2 + C}. (2.49) 

We consider the case in which a particle moves along a constant radius r = ro, but 
precesses around the symmetric axis. The degree of precession is determined by 
the value of C. If ro and C are given, the energy E and the z-component of the 
angular momentum l z are obtained by the two equations, R = and dR/dr = 0. 

Then the energy momentum tensor of a test particle is written as 

T " V = l^JTrf W Hr ~ r ) S (0 - 0(t))^ - *«)). (2.50) 
The source term of the Teukolsky equation is 

/ qauj 

dndt P -*p-\B' 2 + B '* )e -^^t^_0 7 (251) 



where 



B' 2 = - l -p^pL^[p- A L Q {p- 2 -p- l T nn )} 



^p^ 2 L^[p-Yj + {p- 2 T 2 ^ l Trnn)l 



B' 2 * = - ip 8 pA 2 J + [p- 4 J+(p- 2 pT^) 



(2.52) 



with 



p = (r — iacosO) 1 , 
vn 

L s = da H — aw sin # + s cot 6* / 9 r o\ 

sint) l z - JO J 

J+ = d r + — , 

and Q denoting the complex conjugate of Q. In the present case, the tetrad 



* In this section, these constants of motion are those measured in units of [i. That is, if 
expressed in the standard units, E, l z and C in Eq.(2.48) are to be replaced with E/n, l z //j, 
and C/ ii 2 , respectively. 
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components of the energy momentum tensor, T nn , Tmn and T^m, are in the form, 

T nn = ^5(r - r )5(9 - 0(t))6(<p - <p(t)), 
sm9 

T mn = ^5(r - r )S(9 - 9{t))5{ V - cp(t)), (2.54) 
sin (7 

Trnm = ^fS(r - r )6(9 - 0(t))6(<p - <p(t)), 
sin 9 



where 



Cnn = ^- i [E(r* + a 2 )-al z f, 

PP , _2\ _7 H.-_:_ /l/'-TTl ^ 

{-J mm. — 



^^[E(r 2 + a 2 ) - al x ][imne{aE - -4^) + 0(9)}, (2 . 55) 



Cmm = %t[ism9(aE - -%-) + Q(9)] 2 , 
LlA V sin 9' 



2St l V sin^ 
and t = dt/dr. 

Substituting Eq.(2.52) into Eq.(2.51) and performing integration by part, we 
obtain 



00 

r 4 

- £mu> 



2tt 

—00 
2^2 



00 

J ^-^[-^4{p" 4 4(p 3 ^)}Cnnp- 2 p-^(r-ro) 



where 



+ ^j=-{l\s + ia(p - p) sm9S)J+{Crn n p- 2 p- 2 A- 1 5(r - r )} 
V2p 

+ ^Lt{p*S(p 2 p-%}C mn Ap- 2 p- 2 8(r - ro) (2-56) 

- ^p 3 A 2 SJ+{p- 4 J+(pp- 2 C mm S(r - r ))} 

00 



4- ^2 

Ll=d e + auj sin 6> + s cot 0, (2.57) 

suit) 



and 5 denotes _2S^(#(t)) for simplicity. 
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Equation (2.56) can be further simplified by noting that the orbits of our 
interest have the properties, 



6{t + At) = 6(t), (p(t + At) = (p(t) + A(p, 



(2.58) 



where At is the orbital period (of the motion in the ^-direction) and Aip is the 
phase advancement during At. For convenience, we set 



o - 2n o - Av? 



Then the source term reduces to the form, 



T lmu = T one e ik ^ At ~ mA ^ = Q e T one % ~ "»), 



(2.59) 



(2.60) 



where 



uj n =nVt e + mVL^ , (n — 0, ±1, ±2, • • •), 
At 



dte iuJt ~ im ^% mL0 



o 

=A 2 



(2.61) 



(A nn0 + AfnnO + A W mo)S(r - r ) 

+ {(4ni + A mml )5(r - r )} r + {^m2<5(r - r )} rr ] , 



(2.62) 



and the A's are given in Appendix B. Inserting Eq.(2.60) into Eq.(2.5), we obtain 

Zlmcj as 



J lmuj 



n 

Q 



0,7., R« 



Ammo} 



(2.63) 



j pin j2 pin 

a ^mw n r , A x , a n £muj n . 

'l^mnlt/imml; t , 9 %m2 



dr 2 



r=r 



Now, let us consider the orbital integral in Eq.(2.62). To perform them, we 
must know the trajectory of a test particle, i.e., 6{t) and (p(t), but they are not 
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simple analytical functions of t for general trajectories. Thus, we here consider the 
case in which a particle is in the orbit with small inclination. To be specific, we 
introduce a dimensionless parameter y defined by 



y = ^; Q 2 = i 2 z + a 2 (i-E 2 ), 



(2.64) 



and regard it as small. Since Q 2 ~ I 2 and C ~ I 2 + l 2 ] 13] this is equivalent to 
assuming I 2 + I 2 <C I 2 . Note also that we do not need the exact expressions for E 
and l z in terms of tq and C (or y), but only the PNE of them. To the first order 
of y as well as to the P 5 / 2 N order, they are given by 

^ 1 M 3M 2 M 3 / 2 0/1 y, fi , 



Z, =(Mr ) 1 /2 
27M 2 



1 - 



3M y 3MV2a 



(2.65) 



+ 



8r 2 



1 - £) + -(! - 2y) - y) + 0(t; 6 ) 



2r, 



To solve the geodesic equations under the assumption y < 1, we first set 
9 = 7r/2 + y 1 ^ 2 9' and consider the geodesic equation for 9. It then becomes 




(2.66) 



Since the r.h.s of Eq.(2.66) contains only even-functions of y 1 ^ 2 9', we can solve it 
iteratively by expanding 9' as 

0' = 0(0) + !/*(!) + ^(2) + • • • • (2-67) 

[23] 

This is similar in spirit to the method used by Apostolatos et al., who considered 
gravitational waves from a particle in an elliptical orbit around a non-rotating black 
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hole with small eccentricity e C 1. However, here we only consider the lowest order 
solution $(o)- This means we take into account the effect of inclination up to 0(y), 
as seen from the structure of the geodesic equations (2.48). The equation for 9^ 
is 

<dB( 



or dividing it by (dt/dr) 2 , 



f d0 (0)\ 2 _Q A (l n2 , 
f d9 (0)\ 2 _Q 2 n «2, 



(2.68) 



(2.69) 



where 



a = —a(aE — l z ) + 



A(r ) 



{E(r2 + a 2 )-a/4. 



Then the solution is easily obtained as 



(o) = sin(lV); fifl = ^ ^ 



(2.70) 



(2.71) 



where we have chosen 9^ — at t — 0. Thus we have 



7T 



^ = | + l/ 1/2 sin(^). 



(2.72) 



Note that the solution (2.72) implies that the inclination angle 9{ is indeed given 
by 9i = y 1 / 2 in the present approximation. 

Next, we consider the geodesic equation for <p. Taking account of the terms up 
to 0(y), it becomes 



dip k 
dt a 



1 + 



fl z a 2 E 
V k a 

n 



(2.73) 



ftp — y— cos(2f2#t), 



where 



k = -{aE - l z ) + ^^{ E ( r l + « 2 ) - al z }, 
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and 



^ K 1 ^ ( A<zA ^ ^ (l z a 2 E\ . „ , 

% = - + -yn 2 [=^j, n J = _(^— _). (274) 



The solution to Eq.(2.73) with y? = at t = is 



V? = Qpt - y^- sin(2fi^). 



(2.75) 



Note that ^ Qg. This means the precession of a test particle orbit around 
the spin axis of the black hole. Specifically, to the order required for the present 
purpose (see §3 below), we have 



ft, 



M 1 / 2 

3/2 
3/2 



M l l 2 a 3 (M l / 2 a <r \ . 



SM 1 /^ 3a 2 ^. 6 . ' 



(2.76) 



We see that fl v — Hq — > 2Ma/rg for ro — > oo and y — > 0, which is just the 



[24] 



Lense-Thirring precessional frequency. 

Now that we have the solution of the geodesic equations, we can estimate the 
A's in Eq.(2.63). Up to 0(y), they are integrals of the form (see Appendix B), 



At 



W = / ^'-^{qL*. + y 1/2 Qem.J( Q ) + v<& nUn e 



2 

(0) 



(2.77) 

where Q\ mUn for /c = 1 ~ 5 are complicated functions of ro- Substituting Eqs.(2.71), 
(2.75) into Eq.(2.77), and using the approximation, 



(2.78) 
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we find 



Hmuj n 



2tt 

Q f2 2 
+ yjr«-2 - Kl)Q\mu n + (2*n,0 + S n - 2 + S n , 2 )Q 5 emUn 

+ 0(y*/% 

(2.79) 

where 5 n ^ n > is Kronecker's delta. Applying these formulas to As of Eq.(2.63), the 
amplitude Zi mUn is found to have the form, 



J lmu) n 



+ y{Z 2 > 2 5 n , 2 + Z 2 '- 2 5 n ,_ 2 ) + 0{y 3 ' 2 ) 



(2.80) 



where are functions of tq. In principle, algebraic calculations of Z t,J are 
straightforward, but they are almost impossible by hand in practice. Thus to 
avoid trivial mistakes as well as to save time, we have made use of the algebraic 
manipulation program, Mathematica, to obtain Z l,J . We shall not list their explicit 
forms here since they are too complicated and not much insight can be gained from 
them. 



3. The energy and angular momentum fluxes 

In this section, we calculate the energy and angular momentum fluxes to 0(v 5 ) 
beyond the quadrupole formula and to 0(y) in the orbital inclination. From 
Eqs.(2.1), tp4 at r — > oo to the required order takes the form, 

1 2 4 l qaOJ n 

^ = ~ E E E Zern» n ^j^e^*-t) + irn^ (3 . 1} 

It is worth noting that there are symmetry relations, -2S^(9) = -2S^^(n — 9) 
and Z£_ mu _ n = (— l) n Zg mu , n . At infinity, ip4 is related to the two independent 
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modes of gravitational waves h+ and h x as 



^A = \(h + -ih x ). (3.2) 



From Eqs.(3.1) and (3.2), the energy flux averaged over t ^> At is given by 



dE\ x - \ Z £mu 



2 



, 

\ dt I ^ Anuul ^\dt J imn 

£,m,n £,m,n 



In the same way, the angular momentum flux is given by 



2 

dJ z \ _ m \ Z £m^ n _ /dJ z \ _ ^_f^\ (3A) 



(it / f-^ 1 4nuj3 f-^ \ dt J imn f-^ LV n \ dt ) imn 

£,m,n £,m,n £,m,n 

As clear from Eq.(2.80), since we take the square of \Zi muJn \ we only need n — 0, 
±1 modes for the present purpose. Further, since u n = mVL^ + nQg, we only need 
f2<9 at the 0th order of y. The PNE of VL^ and Q,q are given in Eq.(2.76). In order 
to express the post-Newtonian corrections to the luminosity, we define f]£ mn as 

(fL4(f (3 - 5) 

where (dE/dt)^ is the Newtonian quadrupole luminosity: 

d£\ 32/i 2 M 3 32 //ix2 1Q 

^J 7V = ^T = y(m) v • (3 - 6) 

We note that because the mode indices (£, m) here are those associated with the 
spheroidal harmonics, as clear from Eq.(3.1), they do not correspond to the usual 
spherical mode indices. This point should be kept in mind when one attempts to 
interpret the PN corrections to r)£ mn in the language of the standard PN approach. 
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For £ — 2, the results are as follows. If |m + n| > 2 or m + n = 0, i]i mn becomes 
0{v k ) (k > 5). The remaining r\i mn which contribute to the luminosity to 0(v 5 ) 
are given by 

, 107 o o o 4784 4 . 4 428tt «, 4216 . 

%±2 o =1 - + 4tt^ - 6gt; 3 + — w 4 + 2gV - -y-^ 5 + -yjgW 

/ 170 o o o 4784 4 9 4 

+ « -H - 4tcv 6 + 15qv 6 v 4 - lloV 

y \ 21 y 1323 y 

428tt . 13186 = 

H w° qv° 

21 189 y 

1 o 17 , 7T , 17 



/ 1 2 X ' 4 71 5 X ' ^ 

^i=I/(36« +i8 V + li34 9V J' 



1 2 1 3 17 4 1 2 4 ^ ^ 793 k 

772+1 o = — f 9"^ v H o v H f qv 

1 u 36 12 y 504 16 H 18 9072 y 

/ 5 9 1 q 85 4 1 9 4 5n c 13931 r\ . _ N 

+ y I + -qv 3 + v A q V + qv 5 I , (3.7) 

y V 72 8 y 1008 32 y 36 18144 y /' 1 ; 

( 170 2 , o _ o 4784 4 11 2 4 
?72±i±i =j/(l - — ^ 2 + 4™ 3 - 12gt; 3 + — t; 4 + ygV 

428tt . 11078 a 

H of 5 , 

21 189 / 

/ 1 2 1 3 17 4 1 2 4 T 5 745 5 \ 

772 0±1 ^v24 V -12^ ~ 336" + 24 9V + 12 " " 1008 ^ > 

Putting together the above results, we obtain (dE/dt)e = J2mn(dE /dt)i mn for 
£ = 2 as 



r 1277 2 „ o 73 3 / y\ 37915 4 

— = — U tT + 4vn; 3 m; 3 1--)+ v 4 

\dtJ N \ 252 12 y V 2/ 10584 

33 2 4 527 2 4 2561tt . 201575 = / 

+ — g V g + gv 5 1 - - H . 

16 y 96 y y 126 9072 y V 2/1 



(3.8) 



For £ = 3, the non-trivial ?% mn are given by 



1215 2 1215 4 3645tt «, 1215 = 

r?q+q = v v H v qv 

16 3 u 896 112 448 112 H 

( 3645 2 3645 4 10935tt . 3645 «,\ 

+ w tT H w 4 w & H gv^ , 

y V 1792 224 896 112 y /' 

5 4 
%±3ti =^ y, 
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5 4 40 5 / 20 4 100 5 \ 

m±2 = 63 v ~m qv +y {~63 v + m qv )> 

/3645 o 3645 4 10935tt 5 6075 = 

m+2+1 =y f v H f qv 

I6±z±l y\ m2 224 896 224 H 

( 5 2 5 4 5i c 25 r\ 

^i=I/(l6i28« - 3024^+ 8064^ + 18144^)' 

1 2 1 4 71 5 17 ir 

773-1-1 n = v f H f Q"w 

l6±i u 8064 1512 4032 9072 y 

/ 11 2 11 4 Utt = 95 t 

+ y f H f ir H qv" 

y \ 16128 3024 8064 9072 y 

/25 4 80 5 \ 

/ 1 2 1 4 ^ 5 11 5 

7?q n+i =y u v H f qv 

i6 oil y^ 2688 504 -r- 1344 100g y 

The other r]£ mn are of 0(t> 6 ) or higher. Then we obtain 

~"V(l-f)}. (3.10) 



/<Z.E\ /<£E7\ rl367 2 32567 4 16403vr 5 896 
\~dt /3 _ /jv ^1008^ 3024 U + 2016 U ~ ~8T 



For t = 4, we have 



1280 4/ 
?74±4 o = ~^f v (1 - %), 



2560 4 

%±3 T i = Y^i^' ( 3 ' 11 - ) 



»W±2 o = ^ 4 (l-8l/), 
and the others are of 0(t> 6 ) or higher. Hence we obtain 



(dE\ _ fdE\ 8965 4 



V di 74 \dt J N 3969 v 7 



Finally, gathering all the above results, the total energy flux up to 0(v 5 ) is 



24 



found to be 



dE\ fdE\ 



^- 1 



dt I \ dt J 



N 



1247 
"336" 1 



„ s 73 o / y\ 44711 4 

+ 4tto 3 gv 3 1 - v 4 

12 y V 2/ 9072 



33 2 4 527 2 4 8191vr 5 3749 5 



+ 16^ - 96 9 VI/ 



672 



V + 



336 



■gt> 



0-1))- 

(3.13) 

The terms without g agree with those derived in Ref.[15], and the term — 73gt> 3 /12 



also agrees with the previous results. 



[10,12] f 



From Eq.(3.4), the avaraged angular momentum fluxes for t = 2, 3 and 4 are 
calculated to give 



dJ z 
dt 



dJ z 
dt 



dJz 
dt 



dJ z 
dt 

- qv 

— 7TV' 

dJ z 



i-i- 1 ^ (i- y -) + A^ (i-y 

2 252 V 2/ V 2 



N 

61 61 



37915 



12 8 V + 10584 



y 



^(1-^1+ gV I — 



33 229 



,2561 / y 
126 V 2 



qv 



( 22229 27809 



jV 



1367 



1006 



+ 7riT 



dJ, 



,16403 



2016 



(»-!) 



dt 



jV 



8965 
3969 



V 4 (l 



V 1296 
tj \ 32567 
2 



364 



W 4 (l 



g-y 



3024 
88049 9817 



9072 756 



y 



16 32 



y 



(3.14) 



(3.15) 
(3.16) 



where {dJ z /dt)^ is defined to be 



dJ z 
dt 



N 



32/i 2 M 5 / 2 32 / n \2 7 

^77^ = T (&) 



(3.17) 



f Note that Poisson 1121 defines the PN expansion parameter as i>' = (Mf^) 1 / 3 , which is 
related to our v = (M/tq) 1 ^ 2 as v' = v{l — qv 3 /3 + 0(w 6 )) for y = 0. Consequently, his 
Newtonian quadrupole luminosity differs from ours by a factor (w'/u) 10 = (l — 10qw 3 /3 + 
0(v 4 )). This explains the apparent difference between his result — ll<;v 3 /4 and ours. 
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Total flux of the angular momentum is then given by 




dJ z 



dt 



qv 




61 



12 



8191 



672 



We note that the result is proportional to (1 — y/2) in the limit q — > 0. This is 
simply because the orbital plane is slightly tilted from the equitorial plane by an 



From the above results, we find the following features of the gravitational wave 
luminosity: 

1) As argued in §2.1, the quadratic terms in q{= a/M) appear at the v 4 order, and 
the linear terms in q appear at v 3 and v 5 orders. 

2) The coefficients of qv 3 and qv 5 in (dE/dt)2 and (dE/dt)s (and hence in (dE / dt)t t) 
are proportional to l—y/2. Since l — y/2 ~ cos#j, these terms may be regarded as 
proportional to the inner product S • L of the spin angular momentum S and the 
orbital angular momentum L. With this interpretation, our result at the v order 
is consistent with the PN calculation of the spin-orbit terms by Kidder et alf 101 as 

[13] 

well as with the numerical result of perturbative calculations by Shibata. ' 

3) Contrary to the feature 2), the coefficient of the q 2 v 4 term in (dE/dt)2 does 
not seem to be expressible as a simple function of cos 6%. We suspect that a major 
part of it is attributable to the quadrupolar gravitational field around the Kerr 
black hole which modifies the particle orbit. In fact, for y — 0, the 2q 2 v A term of 
7/2±20 in Eq.(3.7) can be explained in terms of the Newtonian quadrupole formula 
as the contribution from the quadrupole moment of the Kerr black hole* However, 
the 1/16 q 2 v 4 term in ?72±io cannot be explained in this way. An inspection of the 

* This was first pointed to us by E. Poisson. 



angle 9{ ~ y 1 ^ 2 , hence dJ z /dt ~ (dJtot/dt) cosflj 
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expanded form of the i = 2 ingoing-wave Teukolsky function given in Eq.(2.35) 
reveals that the q 2 z 2 term at 0(e 2 ) is proportional to m 2 — 4, hence it vanishes for 
m = ±2 while it remains finite for m — ±1. Since this term will contribute to the 
g 2 t> 4 terms in the luminosity, we may interpret the g 2 t> 4 term in r?2±lo as due to 
the curvature scattering in the near zone field. Incidentally, this suggests that the 
coincidence of the coefficient 2 of the g 2 t> 4 term in r/2±20 with the Newtonian calcu- 
lation is rather accidental; naively we would expect the same curvature scattering 
effect to give rise to some additional contribution to the g 2 t> 4 term in ?72±20- In 
any case, our result suggests the existence of a new type of spin-dependent terms 
in the energy flux when a PN analysis beyond the present level is carried out. 



4. Implications to coalescing compact binaries 

In this section, we discuss the effects of PN terms in the luminosity to the 
orbital evolution of inspiraling binaries composed of neutron stars (NSs) and/or 
black holes (BHs). Although our results are valid only in the test particle limit, 
we ignore this fact in the following. 

Since we are interested in the orbits off the equatorial plane, we must consider 
the evolution of C as well as r of the particle as it radiates gravitational waves. 
Although not proved in any sense, let us we assume that the orbit remains quasi- 
circular, i.e., the radius of the orbit is approximately constant for many orbital 
periods. We then would like to see if the inclination angle 9{(= y 1 ^ 2 ) changes in 
time as the orbit shrinks. Since the test particle orbit in this case is characterized 
by two of the four approximate constants of motion, r, C (or y), E and l z , we can 
estimate the change of y with respect to the change of r by equating the energy 
and angular momentum luminosities of the gravitational waves with those lost by 
the particle, 

_ /dE\ dE(r,y) dr + dE(r,y) dy 

\ dt I dr dt dy dt ' ^ ^ 

I dJ z \ _ dl z (r,y) dr dl z (r,y) dy 
\ dt / ~ dr dt dy dt' 
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To the leading order in v and y, this gives 



dlny 61 /M\ 3/2 , A n . 

= - • ( 4 - 2 ) 



dmr 24 \ r / 

Thus y changes only by a small amount during the entire inspiraling stage until r ^ 
10M even for g = 1. Hence the approximation y = const, throughout the evolution 
of the orbit will be good if the orbit remains quasi-circular. Furthermore, by 
adopting a radiation reaction formula which is valid at least in the Newtonian limit, 
it has been numerically found by Shibata' 131 that the evolution of the inclination 
angle is small at least in the low frequency region, r/M ^ 30. This result is 
consistent with the assumption of quasi-circular orbits. Thus we assume y = const. 
in the following. 

Then the total cycle N(ri,rf) of the phase of gravitational waves from an 
inspiraling compact binary from r — Ti (t — U) to r = rj (t — tf) is 

U r f 

where / is the frequency of the wave. Expanding Q^, dE/dt, and dE/dr with 
respect to v — (M/r) 1 / 2 , iV is expressed as 

N 5 M 7 drr^EZoh(q)(M/r) k / 2 EZoCk(q)(M/r) k / 2 

64vr n J M5/2 Er=(A(<7)Wr) fc /2 ' 1 ' ) 

r f 

where the series forms in the denominator and numerator represent the PN cor- 
rections to the Qcp, dE/dt, and dE/dr, that is 

OO Q 

J2b k ( q )(M/r)^ = j^-, 



E^Mfrf'^^, (,5, 

(dE/dt) 



(dE/dr) N ' 
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and the argument q(= a/M) is given to the coefficients b k , c k and d k to emphasize 
the g-dependence of the PN corrections. To the PN order we consider in this 
paper, the PN expansions of and dE/dt are given in Eqs.(2.76) and (3.13), 
respectively. For completeness, here we show the PN expansion of dE/dr: 



dE 
dr 



2^2" 



i 3M ^n ^ 



f + 3^(l 



y) 



r J 



3/2 



My 21 (M\ 

-) -^ q {v) 



5/2' 



(4.6) 



Since the effect of PN corrections in the case of q = has been already- 



studied) 151 here we examine only the effect due to non- vanishing q. For this purpose, 
we introduce the quantity AN^ defined by 



AM n \q) = 



5 M 

64:71 jJL 



drr 3 / 2 ELo b k (q)(M/r) k / 2 ELo c k (q)(M/r) k / 2 



M 5 /2 



ELod k (q)(M/r)^ 



drr 3/2 f Yr=lh{q){Mlr)^ + b n mM/rY^ 
I Tl~l d k (q)(M/r)W + d n (0)(M/r)^ 



n-1 

X 

k=0 



x Ec fc (g)(M/r) fc / 2 + c n (0)(M/r) 



n/2 



(4.7) 



This describes the effect of q corrections at the P n / 2 N order. 

First we consider a NS-NS binary of equal mass M = 1.4M (M/fj, = 4) as a 
typical example. The future laser interferometric gravitational wave detectors such 
as LIGO 111 have good sensitivity in the frequency band between 10 to 1000Hz, so 
we choose T{ = 175M and rj = 8M. A NS of mass 1.4M and radius R ~ 10km 
has q ~ 0.4/P ms , where P ms is the period of rotation in units of a mili-second. 
In this case, the phase N is accumulated at relatively large radii r/M <; 100. 
Hence the convergence of the PN series in the numerators and denominators of 
the integrands in Eq.(4.7) is good enough so that we may expand them further to 
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obtain the approximate formula, 

AN^(q) ^ — - J 
w 64tt fj, J 



5 M } dr (M\ n ~ 3/2 



r f 



M\r ) (4.8) 



x 



{ (b n (q) + Cn(q) - d n (q)) - (6„(0) + c n (0) - d n (0)) } . 



This gives 



AJV< 3 > ~ f ^) (1 - 0.%) , 



2 

in 



ms 



AN (5) „ f*\ (1 _ .3y) . 



Although the y-corrections in the above are valid only for y <^ 1, we expect them to 
be qualitatively valid even for y — 0(1). As mentioned previously, the corrections 
AAK 3 ) and AM 5 ) are due to the spin-orbit coupling, hence replacing y with 2(1 — 
cos^j) in the above formula will give a reasonable estimate in the qualitative sense. 
As for AN( 4 \ although we have not been able to specify the physical meaning 
of it with certainty, at least we may say that the dominant contribution comes 
from the quadrupole moment of the gravitational field around the Kerr black hole, 
as discussed at the end of the previous section. Hence we may also expect the 
replacement y — > 2(1 — cos Q{) to be approximately correct. Thus the inclination 
angle directly affects the values of these phase corrections whenever they become 
important. 

We know three binary pulsars in our Galaxy which will merge within a Hubble 
time, PSR1913+16, [25] PSR2127+11C 1261 and PSR1534+12. [27] Hence these may be 
regarded as a typical target of the gravitational wave detectors. Their rotation 
periods are P ms = 59.0, 30.5 and 37.9, respectively. If we also assume these values 
as typical, we have AiV^ 3 ) > 1 as has been discussed previously) 11 ' 131 while AiV^ 4 ) 
and AN^ are small. However, we know there are several pulsars with P ms ^ 2 
in our Galaxyf 81 for which both AiV^ 4 ) and AiV^ 5 ) exceed unity. Hence it will be 
safe to construct templates which take account of the P 2 N and P 5 / 2 N spin terms. 
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Note however that a main contribution to the correction AA^ 4 ) is due to the 



quadrupole moment of the gravitatinal field, the value of which reflects the pecu- 
liarity of the Kerr black hole. Concerning this point, Bildsten and Cutler' 291 consid- 
ered the quadrupole moment of the gravitational field induced by the quadrupolar 
deformation of a NS due to its rotation and evaluated the phase correction for a 
realistic NS model as 



This result is somewhat larger than our estimate. That is, the effect of the 
quadrupole moment of the gravitational field due to a spinning NS is larger than 
that due to a spinning BH for the same dimensionless spin parameter q. This sug- 
gests that if one body of a compact binary is very rapidly rotating so that we are 



terms which carry information of the orbital parameters such as the spin-orbit 
terms, it will be possible to distinguish a BH from a NS even if the BH has mass 
M B h ~ 1.5M . 

Next, we consider a BH-NS binary composed of 1OM BH and 1.4M Q NS 
[Mj fi = 9.28). Our result has more direct applicability to this case. For simplicity, 
we set Ti = 68M and rj = 6M irrespective of q and y. In this case, the estimation of 
AiV( n ) in terms of the approximate formula (4.8) will not be a good approximation. 
Instead we must use the original formula for AN^ n \ Eq.(4.7), as it is. This is 
because the phase A in the present case is accumulated at smaller r/M than in 
the case of a NS-NS binary, hence the convergence of the PN expansion becomes 
slow. As a consequence it is not possible to derive approximate formulas for AN^ 
as simple as Eq.(4.9). Here we only quote the critical value of the spin parameter q 
above which each correction AN^ (n = 3, 4, 5) exceeds unity. We find AA^ 3 ) 1 
for q ^ 0.01. Hence the correction at this order is important even for a very slowly 
rotating black hole. As for AA^ 4 ) and AN^ , they become larger than unity if 
q 0.2. This indicates that yet higher order PN corrections will be important if 
the BH is rapidly rotating (q ~ 1). 




(4.10). 



able to measure AA^ 4 ) by the matched filter technique, then together with other 
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Summarizing the above analyses, we obtain the following conclusions: 

1) The spin-orbit coupling term at the P 5 / 2 N order is important for the evolution 
of NS-NS binaries with P ms < 2 and BH-NS binaries with q > 0.2. Since the 
rotation of P ms ~ 1 would be the fastest possible period that a NS could have, 
inclusion of the phase corrections up through the P 5 / 2 N order seems to be enough 
for NS-NS binaries. On the other hand, the terms higher than the P 5 / 2 N order are 
likely to be important for BH-NS binaries since the rotation of q ^ 0.2 for a black 
hole seems quite possible. 

2) The q 2 terms at the P 2 N order becomes important for BHs with q <: 0.2 or NSs 
with P ms ^ 2. However, the latter value is based on our formula which is valid 
only for a rotating BH. An estimate base on a realistic NS model gives P ms ^ 5. 1291 
Thus it will be possible to distinguish a small mass BH from a NS if the phase 
corrections to P 2 N order can be detected by matched filtering. 

3) At any order of PN corrections, the effect of a finite inclination angle to the 
number of the phase cycles must be taken into account whenever the spin terms 
become important. 

5. Summary 

In this paper, we have performed a post-Newtonian calculation of the gravi- 
tational waves from a particle of mass /i orbiting around a rotating black hole of 
mass M (/i <C M). We have considered the orbit of a constant coordinate radius 
r = tq but with small inclination angle 6{ ~ y 1 ^ 2 , where y is a non-dimensional 
parameter proportional to the Carter constant of the orbit. 

We have formulated the post-Newtonian expansion of the Teukolsky equation 
and its source term in terms of a small expansion parameter v = [Mj tq) 1 / 2 accurate 
up through 0(v 5 ) (P 5 / 2 N order). We have not directly dealt with the Teukolsky 
equation but first formulated a method to obtain the homogeneous solution to 
the Sasaki-Nakamura equation, which is a generalization of the Regge- Wheeler 
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equation for the Schwarzschild black hole, by expanding it in powers of e = 2Moo, 
where u is the frequency of a gravitational wave. In particular, to 0(e), we have 
obtained the ingoing-wave radial functions for arbitrary spherical indices (£, m) in 
closed analytical form. Then we have obtained all the necessary radial functions 
to the required accuracy and transformed them to the corresponding Teukolsky 
radial functions, which have been used to construct the Green function. Further, 
we have formulated the post-Newtonian expansion of the source term for circular 
orbits with small inclination angle. Assuming ?/ < 1, we have analytically solved 
the geodesies of a particle accurate to 0(y) and obtained the source term to the 
required accuracy. We have used these results to integrate the Teukolsky equation 
and derived the formulas for the gravitational energy and angular momentum fluxes 
which are correct to 0(t> 5 ) and to 0(y). 

Based on thus obtained luminosity formula, we have estimated the accumu- 
lated phase N of gravitational waves from inspiraling binaries, assuming the orbit 
remains quasi-circular. Specifically we have considered a NS-NS binary of equal 
mass 1.4M and a BH-NS binary of masses 1OM and 1.4M , which will be typ- 
ical targets of the near-future laser interferometric gravitational wave detectors. 
We have found that if the rotation of a neutron star is moderate, say P ^ 20 ms, 
only the phase correction at its leading P 3 / 2 N order will be important. However 
we have also found that if one body of a binary is a rapidly rotating NS (P ^ 2 ms) 
or a rotating BH (q = Jbh/M 2 > 0.2), the phase correction of AN > 1 will be 
caused by the spin terms at P 2 N and P 5 / 2 N orders. Furthermore if one body is a 
rapidly rotating BH (q ~ 1), it is expected that the higher order corrections such 
as P 3 N ones become important. In all these cases, when the phase correction at a 
certain PN order become significant, that due to a non-vanishing inclination angle 
at the same PN order becomes equally important. 

The above conclusions imply that it is desirable to evaluate yet higher order PN 
spin corrections to the gravitational wave luminosity. As for the inclination of the 
orbit, since we expect the expansion in powers of y to be valid for y ^ 1 (0j ^ tt/4), 
it will be meaningful and useful to calculate the higher order corrections in y along 
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with higher order PN calculations. These problems are left for future work. 
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APPENDIX A 

In this Appendix we show the potential functions F and U of the SN equation 
(2.10). Details of the derivation are given in Ref.[18]. 

The function F(r) is given by 

F{r) = ^-At* (AT) 

where 

r] = c + ci/r + c 2 /r 2 + c 3 /r 3 + c 4 /r 4 , (A.2) 

with 

c = -YliujM + A(A + 2) - \2auj{auj - m), 
ci = 8ia[3au — \(au — m)}, 

c 2 = -2AiaM(auj - m) + 12a 2 [l - 2(au - m)\ (A.3) 
c 3 = 24ia 3 (auj - m) - 24Ma 2 , 
c 4 = 12a 4 . 

The function U (r) is given by 



U ( r ) = , 2 A \ 2 + ° 2 + 4^2 ~ FG ' ( A - 4 ) 
(r 2 + a 2 ) 2 r 2 + a 2 



34 



where 



G _ _2(r-M) rA 

r 2 + a? (r 2 + a 2 ) 2 ' 



K(3 T ^ , 6A 
a = -1-^2" + 3 ^ + A + 72" > 

/3 = 2A (-iff + r — M — ^) . 



APPENDIX B 



In this Appendix we show the A's in Eq.(2.60). 



ImmO 



A- 



^mm2 



At 



A„ n0 = ^ J dte^-^Cnnp-'p-'L+ip-'L+ip'S)}, 




At 



ArnnO = / dte^^ C WnP ^ [(L+ S) (^+P + p) 



.K 



a sin 6 S—(p — p) 



At 



At 



A-mn 1 — 



%- [ dte^ t - im ^p- 3 Crnn[LtS + tasm6(p-p)S], 
V^A J 



At 



mm 1 — 



J dte wt ~ im ^p- 3 pC m mS 



K 



At 



-±= f dte^- im ^p- 3 pC m mS, 
V 2-k J 



where S denotes -2<S^. 
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APPENDIX C 



In this Appendix, we show that the asymptotic form of Xff^ at 6 C 2 € 1 
has the form, 

X tL = z l+1 \0(l) + eOiz- 1 ) + e 2 0(z- 2 ) + e 3 0(^ 3 ) + . . . . (C.l) 



From Eqs.(2.10) and (2.18), the equation for ^ m becomes 
'Ad 2 „ d 



r 



+ C1— + C2 



2 dr 2 dr 



itra = 0; 



C 2 = ^ ( -C/i + 2 - ^ - A0 2 r - A, r Fi + iA0, r Fi - 2iMcu ] , 



where F\ = r) jr /r), and 0, 77 and C/i are given in Eqs.(2.17), (A. 2) and (A. 5), 
respectively. 

Let us now examine the behaviors of the coefficients C\ and C2. By dimensional 
consideration, they must be of the forms, 

C k = ^f k (2M/r,ujr) (A; = 1,2), (C.3) 

where fk(y,z) are dimensionless functions of their arguments. We then easily see 
that fk(y, z) are regular at y — 0, since both C\r and C2?" 2 have well-defined limits 
as M — > 0. Furhtermore, by examining the behaviors of Ci and C2 as — ■> 0, we 
find they are also regular in this limit. Hence -2) are regular at z = as well. 
Thus may be expanded as 

00 

/*(v,*) = £/£ n) (*)v B , (C4) 

n=0 
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in) 

where (z) are regular at z — 0, hence may be futher expanded as 



4 n) w = £ 4"""^ 



(n,m) r 



(C.5) 



m=0 



Note that ff\z) = 2 and ff\z) = -£(£ + 1) + -z 2 , which are the coefficients 
appearing in the lowest order differential operator L^ ) given by Eq.(2.20). 

Taking the above consideration into account, scaling r to z = oor in Eq.(C2), 
and noting that y = 2M/r = e/z, we have 



\ z 4 z 2 J dz 2 



vn=l 



1 d 



z dz 



?(n), 



\n=l 



6' 

7 r 



ilm — 0- 



(C.6) 



Accordingly, if we expand £ £m as & TO (e; z) = J2™=0 ^^imi*) and set ^fl( z ) = 3l ~ 
2^, we can easily show recursively that the asymptotic behavior of £^ at z 1 is 



t (n) 



0(z- n+e ). 



(C.7) 



Then the conversion of £f m to A 7 "^^ by Eq.(2.18) yields the result Eq.(C.l). Since 
e = 0(f 3 ) and z — 0(v), in order to calculate the energy and angular momentum 
fluxes up to 0(v 5 ) beyond Newtonian, we conclude that the necessary power series 
formulas for around z = are those of n < 2 for £ = 2, n < 1 for £ = 3 and 

n = for £ = 4. 
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APPENDIX D 



In this Appendix we describe the expansion of the spheroidal harmonics -2-S^ 
and their eigenvalues A in powers of au. Since we are interested in the energy and 
angular momentum fluxes to 0(v 5 ) in this paper, we need the expansions of -2<S^ 
to 0{auo) for both £ = 2 and 3, but those of A to 0{{auo) 2 ) for £ = 2 and to O(aw) 
for I = 3, while we only need the lowest order formulas for £ = 4. 

The spheroidal harmonics of spin weight s = — 2 obey the equation, 



1 dr. 



sin9 d9 



2 2-2/1 

-a uj sin 



(m - 2cos6>)' 



sin 2 # 

+ 4aa> cos 9 — 2 + 2mauo + A 



(D.l) 



-2^ m = U. 



We assume they are normalized according to Eq.(2.2). We expand -iS^ and A as 



-2SZ = - 2 Ptm + awS® + {au) 2 S^ + 0((a^) 3 ), 
A = A + ac^Ai + a 2 cj 2 A 2 + 0((acj) 3 ), 



(D.2) 



where _ 2 P£ m are the spherical harmonics of spin weight s = —2, Ao = (£ — l)(£ + 2) 
and Ai = -2m{£ 2 + £ + A)/{£ 2 + £) m 

Inserting Eq.(D.2) into Eq.(D.l) and collecting the terms of equal orders in auo 
and (aco) 2 , we obtain 



£o S tm + X o S im = "(4cos^ + 2m + Ai) ^Plm , (D.3) 
C S^ + \ S^ = -(Aco S 9 + 2m + \ 1 )S^-(\ 2 - sin 2 9) _ 2 P lm , (D.4) 



where £q is the operator for the spin-weighted spherical harmonics, 



£d[-2-Pfo 



I d r . „ d \ (m-2cos#) 2 



sin 9 d9 

-A -2-Pfon 



sin6> 



-) 



sin 2 # 



(D.5) 
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First we solve Eq.(D.3) for sj^. Setting 

c(i) \^ a n 



(D.6) 



we insert it into Eq.(D.3), multiply it by - 2 P(i m an d integrate it over 9. Then 
noting the normalization of the spheroidal harmonics, we have 



c ir, 



{£' -!){£' + 2) -{£-!){£ + 2) 



0, 



J- 2 P £ , m cos9- 2 P £md cos9, 

£' = £. 
(D.7) 



Hence cj is non-zero only for £' — £ ± 1, and we obtain 



J+l _ z 

tm ( o j 1\2 



(£ + 3)(£ - !)(£ + m + 1)(£ - m + 1)] 1/2 



(2£ + l)(2£ + 3) 



2 

'£2 



(l + 2)(l-2)(l + ro)(l-m) 
(2£+ 1)(2£- 1) 



1/2 



(D.8) 



Next we consider A2. We set = J2 e , d\ m -iPi'm and insert it into Eq.(D.4). 
This time we multiply it by - 2 P( m an d integrate it over 9, noting that d\ m = 0. 
The result is 



A< 



= -4 J - 2 P£m cos 9 sf^d cos 9 + J - 2 P£ m sin 2 9 - 2 Pe m d cos 9 

= -2(£ + l)^* 1 ) 2 + 2£(cJ- 1 ) 2 + 1 - J - 2 P tm cos 2 - 2 P lm dcos9. 



(D.9) 



The last integral becomes 

/" D 2 fl D , fl 1 , 2(£ + A)(£-3)(£ 2 + £-3m 2 ) 

/ -2-P/m COS 9 -oPfmd COS 9 = - + -- ; - V — L 

J 2 lm 2 lm 3^3 1)(2£ + 3)(2£ - 1) 

We need A2 only for £ = 2. In this case, the final answer becomes 

90 - 10m 2 



(D.10) 



A 2 (^ = 2) 



189 



(D.ll) 
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